Comments on an alternative theory for the accelerating universe 
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We utilize the recently released supernova data of Riess et al. on the acceleration of the universe to 
determine how the parameters of a particular modification of gravity, the DGP 5-dimensional theory, 
are changed when compared to determinations using earlier data sets. We show the parameters of the 
theory are now very tightly constrained and at very reasonable values. For example, the curvature 
parameter, Qk, in contrast to an earlier analysis which used a smaller data set, is now consistent 
with zero, and the redshift at which the universe starts to accelerate is near unity. 

PACS numbers: 98.80 -k, 98.80.Es 



1. INTRODUCTION 

Since the discovery that the expansion of universe is accelerating [l|, most of the work has been in the 
direction of describing this acceleration by a negative pressure component - the dark energy - within the usual 
theory of gravity. Alternative theories try to eliminate the dark energy by modifying the theory of gravity itself 
either by adding extra dimensions in which gravity, but only gravity, propagates [3| or by modifying Einstein's 
equations in four dimensions Q • Here we will consider one such alternative theory, the so-called DGP model 
where gravity is effectively 4-dimensional at short or moderate distances but 5-dimensional at large distances. 
An extensive review of the DGP theory has been given by Lue [f|. The supernova data 0] have been fit to 
this model with the result that the curvature was required to be very large [1, Q . As shown by Maartens and 
Majerotto [Ttj, once the comoving distance to the surface of last scattering given by WMAP [111 . Il2| | is included 
as a constraint on the data, the matter density and curvature agree with other measurements such as the latest 
results from WMAP [l2j]. Here, we examine the implications of using the latest supernova data set of Riess et 
al. [HI to fit the parameters of the DGP model. 

It is known that the form of the Hubble parameter in the DGP theory is indistinguishable from a dark energy 
model H, We check that the equation of state of this effective dark energy, w(z), is well behaved - it varies 
from —0.8 at small redshifts to —0.5 at large redshifts. It has been argued, however, that the DGP theory allows 
values of the effective w(z) that are less that -1 [HEl]; we will discuss this point briefly. 



2. FORMALISM AND RESULTS 



The DGP theory has only one parameter, r c , the crossover radius where the theory changes between a region 
that is effectively 4-dimensional to one that is fully 5-dimensional, defined by 
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where M5 is the 5-d Planck mass. If we introduce a density parameter associated with r c as 

1 



fi r = 

then the Hubble parameter, H(z), is given by [1, 03] 
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where 2 is the redshift and Hm and fifc are the matter and curvature density parameters. There are two possible 
phases denoted by e = ±1. For obvious reasons we will use the self-accelerating e = +1 phase, but we will have 
occasion below to refer to the other phase. From ([3]) we must have 
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We will treat Q r and m as the independent parameters. 

To confront this theory with the data we need the luminosity distances which are given by the usual integral of 
1/H. Prior to the recent supernova data release of Riess et al. [13j |. the largest data sets available were the 164 
supernova events collected in Ref. and Ref. fl7| . th e sum of which we call the essence data, and the legacy data 
set of Ref. (l8| . Fits of the data in Refs. and to DGP model were performed by Maartens and Majerotto 

1C I . We reproduce their results (including the data in Ref. [TtJ ) for comparison with the fits to the data in Ref. 

l3. 



For the essence data, a x 2 fit varying O r , Oj,; and an offset gives 
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where the errors are one sigma. This agrees with the results of Ref. Q which used the gold data and looks 
quite reasonable until we use Eq.Q to find fl k . The result is an unattractively large value 
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Alternately use of the legacy data of Ref. [18j gives somewhat better numbers but with larger errors, 
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However, there is one more luminosity constraint that we can add to the data sets and that is given by the 
comoving distance to the surface of last scattering measured by WMAP. In Ref.[ll| this is called the ' ang ular 
size distance to the decoupling surface', cLa, and it is what Weinberg calls the proper-motion distance [l9|]. In 
particular, we include this data point as a (Gaussian) prior on the parameter (2p| . l2l| Rcbm = v^Ta/^a = 
1.70 ± 0.03 at z = 1089 [Hj. Including this one additional point with its very long lever arm changes the results 
of Eqs.©, ©, and © to 
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The other data set which gave Eqs. 



, and (|10[) now gives 
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FIG. 1: The 68.3%, 95.4%, and 99.7% confidence con- 
tours in the Qm — Q r plane are shown for the Hubble 
parameter of Eq. <[3j) using the essence data without 
the CMB point. The cross denotes the \ 2 minimum 
located at Qm — 0.348, Q r = 0.254, and the line is 
where Q k = 0.0. 



FIG. 2: The 68.3%, 95.4%, and 99.7% confidence con- 
tours in the Qm — Q r plane are shown for the Hub- 
ble parameter of Eq. © using the R06 data with no 
priors. The cross denotes the \ 2 minimum located 
at Qm = 0.358, Q r = 0.237, and the line is where 
Q k = 0.0. 



These results are much better. The problem of the large negative Q k has gone away and the errors are much 
smaller. The Qm value for the essence data looks a bit small but is consistent with the most recent WMAP 
results. In the case of the legacy data, the value of Qm is low. 

Repetition of the these calculations with and without the CMB point using the 206 supernovae of Ref . [HI , 
which we will refer to as the R06 data, gives 
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and 
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respectively. When performing these fits, we followed the procedure employed in Ref. [l3j and did not include 
the 24 supernovae with z < 0.024. Again, the R06 data without the CMB point favors a large value of Qm and 
a large negative value of Q k - When the CMB point is included, Qm value is quite reasonable and Q k is small. 
The confidence contours for the essence data and the R06 data without the CMB point are compared in are 
FigsQ]and|H anc ^ ^he corresponding comparison with the CMB point included is shown in Figs. [3] and 0J What 
is particularly striking is that in either case the inclusion of the CMB point tightly constrains the parameters 
Qm and Q r along the Q k = line, with the R06 data preferring a wider range of values for Qm- 

In Ref. [l(| , fits to the DGP model using the gold data of Ref. Q included a prior on the baryon acoustic 
oscillation parameter (BAO) [H, HI A = 0.469 ± 0.017 as well as the CMB constraint and, in the Ref. [13 ] 
analysis of the ACDM model, fits combining priors on Qm = 0.28 ±0.03 and A (weak prior) and A and the CMB 
point (strong prior) were discussed. The implications of the weak and strong priors on the DGP model using the 
R06 data are tabulated in Tableland illustrated in Figs. [5] and [6l In the case of the weak prior, the inclusion of 
a prior on mass parameter has the effect of orienting the confidence contours in such a way as to favor a narrow 
range for values for Qm as well as significantly reducing the size of the allowed parameter space. Interestingly, the 
strong prior, which only involves the baryon acoustic oscillation and CMB constraints, shifts the most probable 
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FIG. 3: The 68.3%, 95.4%, and 99.7% confidence contours 
in the Qm — Q r plane are shown for the Hubble parameter 
of Eq. using the essence data and the CMB prior. The 
cross denotes the \ 2 minimum located at Qm = 0.191, Q r = 
0.154, and the line is where Hj. = 0.0. 



FIG. 4: The 68.3%, 95.4%, and 99.7% confidence contours 
in the Qm — plane are shown for the Hubble parameter of 
Eq. 10 using the R06 data and the CMB prior. The cross 
denotes the \ 2 minimum located at S)jf = 0.231, Q r = 
0.142, and the line is where Q^ = 0.0. 
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TABLE I: The results for fits to the R06 data to the DGP model are shown for various combinations of priors. MAT denotes the 
matter density prior Qm = 0.28 ± 0.03, BAO denotes the baryon acoustic oscillation prior A = 0.469 ± 0.017 and CMB denotes the 
WMAP prior R C mb = 1-70 ± 0.03. Note that for all cases x 2 /DOF < 0.91. 



value of S1m to about the same value as obtained with the weak prior, in addition to orienting the confidence 
contours in a narrow range along the Sl^ = line. Including the extra 24 R06 events with z < 0.024 changes the 
f2 r , fJj\f, Slfe values in Table|T]by at most 0.004 (and usually much less). The \ 2 increases by about 13 in each 
case. 

Given these very precise values for the parameters what do they predict for the redshift where the universe 
starts accelerating? The deceleration parameter can be written 



-! + (! + 2-) 
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where the prime indicates differentiation with respect to z. Using Eq. ([3]) this gives, after some amount of work, 

?(*) = ^ (16) 

where, with s = Q r /QM, 
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FIG. 5: The 68.3%, 95.4%, and 99.7% confidence contours 
in the Slj\j — $l r plane arc shown for the Hubble parameter 
of Eq. lO using the R06 data and the and the 'weak' prior 
consisting of the the priors on Qm and the baryon acoustic 
oscillation parameter A. The cross denotes the \ 2 minimum 
located at Qm = 0.281, Q r = 0.193, and the line is where 

n k = o.o. 




FIG. 6: The 68.3%, 95.4%, and 99.7% confidence contours 
in the Qm~ plane arc shown for the Hubble parameter of 
Eq. 10 using the R06 data and the 'strong' prior consisting 
of the priors on the baryon acoustic oscillation parameter A 
and the CMB parameter Rqmb- The cross denotes the \ 2 
minimum located at f2jv# = 0.284, f2 r = 0.125, and the line 
is where Qfc = 0.0. 



Since Qk > the denominator is always positive. Then it is not too hard to see that the redshift where N is 
zero, z* , is given by 

1 + z* = 2s? (19) 

For the essence parameters in Eqs. pT]) we find z* = 0.862 ± 0.114. Using the legacy parameters of Eqs. (fT2")l the 
value is slightly higher, z* = 1.082 ± 0.137. The result for the R06 parameters, Eqs. (H]), is z * = °- 701 ± 0- 100 - 
Note that Eq. (|19|) also tells us that s must be larger than | for acceleration to have started by today. 

Now since the Hubble parameter of this theory is known to be identical with that of dark energy we want to 
see what form the dark energy equation of state takes. We will compare the dark energy with the general form 
of Eq.Q which allows for both phases. For a dark energy description, the Hubble parameter is given by 



H\z) 



fi fc (l + z) z + M (1 + zY + n x f{z) (20) 
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where VLxf{z) is the dark energy density and f(z) is given in terms of the dark energy equation of state by 



f(z) = expjs^^a + ^x)) 



Comparing Eqs.© and (|2U)) we sec that 



n x f(z) = 2Q r ± 2ey/Or^/n r + fl M {l + zf (22) 
From Eq. (f2"Tj) we see that w{z) is given by 

M*) = -1 + \( 1 + z )j^j ( 23 ) 

where, as above, the prime indicates differentiation with respect to z. Applying this equation to Eq. (|22p we get 
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where s was defined above. For both phases the second term on the righthand side of Eq. (f2~4"|) is positive for all 
s > so w(z) is always greater than —1 . With e = +1 both of the above fits give uu(z) m —0.85 at z = and 
w(z) — > —0.5 at large z. 

It has been claimed that the phase with the minus sign in (jHJ) can give w(z) less than —1 [Tsl. [l6j. Rather than 
Eq.© the authors of Ref.pjJ use 



= -VO^+ v ( 1 + V^-) 2 + ^M(l + z) 3 -n M (25) 



Hjz) 
H 

Note that the righthand side of Eq. (|2"5|) is unity when z = 1 so O r and S]« can be treated as independent 
variables. Now we can compare Eq. ([25)1 with Eq. (|2"0]) with S\. set equal to zero, determine fix f(z), and use 
Eq. ([23|) to find w{z). The general expression is not very pretty but at z = it reduces to 

which is less than —1 for Q,m < 1. We tried to fit the data to this form of the Hubble parameter although we 
didn't expect much success since Eq. (|25p does not have the self-accelerating phase. We left out the WMAP 
constraint because Eq. ([23]) is only intended to work for small z. The data tried to force f2 r and £Im to be large, 
il r > 10 and Q,m > 1, in which case w(0) is not less than —1. When we added a prior of Qm = 0.30 ± 0.03 then 
Q r is very small, consistent with zero, 0.015 ± 0.023, but w(0) is indeed less than —1. 



3. CONCLUSIONS 



We have considered the DGP theory of gravity and shown that if the expression for the Hubble parameter given 
by Eq. ([3]) holds for all redshifts back to the surface of last scattering so that the comoving distance measured 
by WMAP can be added as a constraint on the supernova data, then the parameters Q r , CIm and Qk can be 
precisely predicted. Our results are given in Eqs. (fTTj) for the essence data, Eqs. (fl2|) for the legacy data, and in 
TablcQ] for the R06 data. The effect of adding the CMB point is illustrated by the confidence contours for the 
essence and R06 data shown in Figs. [T]-@]and the effect of using the weak and strong priors with the R06 data 
is shown in Figs.[5]and[ni 

This theory is consistent with a flat universe which began accelerating at a redshift equal to, or slightly less 
than, unity as given by Eq. and is identical, as far as the deliberations of this paper are concerned, with a 
model of dark energy whose equation of state is given by Eq. (|24|) with s w 0.6 — 1.1. 
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